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Abstract 

This research is an attempt to reduce the gain variation, 
noise figure and to improve the gain broadness of EDFA by 
modified mathematical modeling of EDFA for 96 DWDM 
systems. An improved simulated model of EDFA is specially 
designed after all the major impairments are taken into 
account like noise, ASE fiber length, input pump and signal 
power. The mathematical model of EDFA has been proposed 
by improving the rate equations of EDFA. Variation of gain 
versus wavelength has been analyzed with and without ASE. 
This research claims to support 96 DWDM channels at a 
channel spacing of 0.8 nm, with a gain of 23.8 dB, ASE of 0.9 
dBm for 6 m EDF length for 1479 nm-1555 nm EDFA. 
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Introduction 

The main parameters in the design of an EDFA 
include the fiber glass material, the waveguide 
characteristics of the fiber, the erbium concentration 
profile, fiber length, pump sources and any active or 
passive components such as couplers, isolators, filters 
etc. The design of the amplifier depends on the 
intended application, however, the primary design 
goals are high gainand output power, low noise figure, 
the flat gain spectrum, reliability, etc. The importance 
of each of these parameters depends on the particular 
application type. For example, a single-channel inline 
amplifier requires high power and low noise figure, 
but the gain flatness is insignificant. On the other side, 
in the design of a DWDM amplifier, broad gain profile 
is the most critical parameter. Power amplifiers 
require high gain and output power and may not have 
constraint on the noise figure and gain flatness. 
Analog amplifiers are the most demanding since all 
the described parameters are important to amplify the 



signal keeping the waveform as much undistorted as 
possible. So, one of the main issues when EDFA 
amplifies multi-wavelength signals is the gain 
broadness in the used wavelength range. Previous 
works showed that maximum gain flatness is achieved 
when average inversion population is in the range 
from 0.75 to 0.8, which implies that when EDFA 
operates in the saturation regime, it is not suitable to 
work in DWDM system. Various gain broadening 
techniques of EDFA are available. The aim of this 
work is to optimize EDFA parameters in order to 
obtain high gain broadness and low noise figure. 

Modeling of EDFA for Dynamic Effects 

The basic operation of EDFA is shown in Fig. 1. 
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FIG 1 THE BASIC OPERATION OF EDFA 

The three level energy diagram of EDFA is considered 
in which the signal gain is achieved by a metastable 
population of excited ions by emission decay from a 
higher pumped state. In this section, modeling of 
EDFA is proposed using improved rate equations of a 
three state EDFA by considering forward ASE. Figure 
(2) describes three level energy diagram of EDFA with 
various energy transitions. The three population states 
of Er+3 are ground state (g) with population density of 
ng, the metastable state (m)with population density of 
nm which is related to signal frequency and excited 
state (e) population density of ne which is related to 
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pump frequency. The total population of erbium ions 
at these three states is assosiated with the population 
density (q) of erbium ions. 

n e +n m +n g = p (1) 
The important transition in this case is from ground 
state to metastable state, the energy difference among 
these states corresponds to about 1550 nm. Let Pge be 
the pumping rate from ground state (level g) to excited 
state (level e), Peg and Pern be the stimulated 
emission rate from excited state to ground state and 
metastable state, respectively. It is assumed that Peg is 
not considered as an important transition. There are 
two types of transitions occuring in excited state, one 
of which is radiative transition and the other is non- 
radiative transition. The radiative transition from 
excited state is further of two types i.e. upto 
metastable state (level m) and upto ground state i.e. 
Aem(r) and i.e. Aeg(r) respectively. It is also 
considered that the transition is mainly non-radiative, 
which implies that non-radiative transition (Aem(nr)) 
radiative transition(Aem(r)Aeg(r)). Let the rate of 
stimulated absorption and emission be Sgm and Smg 
respectively. The rates of spontaneous emission from 
metastable state are also radiative and non-radiative in 
nature, at this level radiative transition (Amg(r)) non- 
radiative transition (Amg(nr)). The non-radiative 
transition from excited state and radiative transition 
from metastable state are considered as ne/x' and nm/x, 
separately, where x' and x are the respective transition 
rates. The energy transistion of EDFA is shown in Fig. 
2. 
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FIG 2 ENERGY TRANSITIONS IN A THREE LEVEL EDFA 

Therefore, the improved rate equations of the three 
states for this proposed model are shown in Eq. 2-4: 
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Considering the assumptions on the probability of 
occurrence of transition from excited state to 



metastable state (which is more as compared either to 
that from excited to ground state or from metastable 
state to ground state) and total population inversion, 
the equation (4) becomes 

§ t = (^m n g — Smg n m ~ + Pge n g ) (5) 

The equation (5) can be represented by means of block 
diagram shown in fig.3. 



D.l.T 



Sni >0 



A 
D 
D 
E 
R 



FIG 3 BLOCK DIAGRAM REPRESENTING POPULATION 
DENSITY OF ER+3 IONS IN METASTABLE STATE 

The signal and pump rates(Sgm, Smg) which are 
dependent on cross-sectional are of core (A), 
absorption and emission cross section data (a) of 
erbium fiber and confinement factors for pump and 
signal waves (T) can be defined as: 



£ r s <ig m p s 

^gm - J 
q r s a m g P s 



(6) 
(7) 



The pumping rate from ground state to the excited 
state can be represented in terms of pump power, 
confinement factor for pump by equation (8) 

_ r n (Tap Pr) 



(8) 



On considering scattering losses of the erbium doped 
fiber negligible, the rate equation for the population (5) 
can be rewrite as: 



dJ t=P s (0,t) 



P s (L,t) + P P (0,t)-P P (L,t)- N -f{9) 



Equation (9) represents the population density Nm as 
a function of position along the fiber without 
considering the ASE. The last term in R.H.S. of 
equation (9) is the decay in metastable state 
originating from spontaneous emission, and the rest of 
the terms combines both absorption and emission due 
to fields. Regarding the effect of co-propagating ASE, 
the equation (9) becomes: 

3Nrn = P s (0,t) - P s (L,t) + P P (0,t) 



dt 
-Pp(L,i) 



■ + ^e(0,0 



P; sE (L,t) (10) 
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Where, ASE power, pump power and Signal Power 
are shown in Equation (11) 



= 2hvAv 



(G-1)+-^JV„ 
pA 



Where, B P = 



P P (0,t) 
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(11) 



The equation (11) is considered according to Bononi 
and Rusch's notation. The improved rate equation of 
metastable state is shown in equation(12): 

dJ ^= P s (0, t) [1 - exp(B s N m - C s )] + P P (0, t) [1 - 
exp(£ P W m - C P )] - N -f- P; SE (0, t) [l - 2hvAv(G - 1) + 
2hv^AvN m ] (12) 

The equation (12), the primary equation to describe the 
dynamic gain effects in the three level EDFA including 
ASE, can be modeled using SIMULINK as already 
proposed in. In this work, the improved non-linear 
ordinary differential equation is implemented as a 
MATLAB Simulink model(shown in figure (4)), 
including the impact of ASE. Algorithm and Flowchart 
used for simulation of proposed SIMULINK model of 
EDFA is as follows: 

Algorithm I Algorithm_Sim_EDFA 



STEP I: Initialize ng, nm and ne (Er+3 ion densities at 
ground, metastable and excited states), A(area), 
L (length of Fiber), Pp&Ps (Pump and Signal Power), 
Ap&As (Pump and Signal Wavelength), ASE power 

STEP II: ng, nm, ne and Invariable 

STEP III: Calculate ion of Er+3 in metastable state 
and length of EDF 

STEP IV: Calculate Gain, ASE Power and Output 
Power with respect to Length of Fiber 

STEP V: Calculate Optimum Length of fiber for 
maximum gain and minimum ASE power 

STEP VI: Plot gain, NF for optimum length of EDF 
w.r.t.wavength 

STEP VII: Goal Achieved 



By means of Algorithm I and figure (4), the graph is 
plotted between gain and wavelength (without 
considering ASE) for EDF length ranging from 4 m to 
14 m as shown in figure (5)which implies that the 
length 14 m can be chosen as optimum length, because 
at this length peak gain of 32 dB is obtained while for 6 
m its value is 26.5 dB and for 4 m its value is 20 dB. 



The main goal of this work is to define the optimum 
value of EDF length so that gain broadness is 
maximum in the largest possible wavelength range. 
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FIG 4 IMPROVED SIMULINK MODEL OF EDFA 
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FIG 5 GAIN VERSUS WAVELENGTH FOR DIFFERENT FIBER 
LENGTHS WITHOUT ASE 
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FIG 6 GAIN VERSUS WAVELENGTH FOR DIFFERENT FIBER 
LENGTHS WITH ASE 

To obtain the optimum fiber length, gain of EDFA is 
plotted as shown in figure (6) by taking the effect of 
co-propagating ASE into consideration for length of 
EDF ranging from 4 m to 14 m. When ASE is included, 
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the optimum length for gain is determined by 
maximizing the signal gain which is a function of the 
length. It has been observed that the gain peak is 
reduced from 32 dB to 26 dB for 14 m length, from 26.5 
dB to 23.8 dB for 6 m length and from 21 dB to 15.5 dB 
for 4 m length. It has been clear from figure 7 that 
there is a very insignificant difference in gain peaks for 
14 m length and 6 m length. So, it is concluded from 
figure 6&7 and Algorithm I that the fiber length of 6 m 
is considered as the optimum length for pump power 
of 17.6 dBm. 

The simulation results presented in figure (5) and 
Figure (6) can be summarized in Table 1. 

TABLE 1 SIMULATION RESULTS 
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23.3 
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12 


30.1 


12 


25.3 


9.1 


14 


3 2 


14 


26.5 


11.8 



The results show that EDF length of 6 m can be 
considered as optimum length and for this length the 
gain and noise figure are plotted with respect to 
operating wavelength as figure(7) which shows 
improvement in broadness of gain spectrum. 
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FIG 7 GAIN AND NOISE FIGURE VERSUS WAVELENGTH 

Conclusions 

The improved SIMULINK model enhances the 
broadness of gain spectrum of EDFA. The optimized 
values of population densities, cross-section areas of 
absorption and emission spectra, length of EDF, pump 
powers are considered. Though ASE is considered, 



even then broadness of gain spectrum gets improved, 
because a compromise for optimal length at which 
ASE is negligible is made, the value of maximum gain 
is lowered. 
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